Equilibrium Charge Distribution on Annealed Polyelectrolytes 
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Monte Carlo simulations are used to study the non-uniform equilibrium charge distribution along a 
single annealed polyelectrolyte chain under 9 solvent conditions and with added salt. Within a range 
of the order of the Debye length charge accumulates at chain ends while a slight charge depletion 
appears in the central part of the chain. The simulation results are compared with theoretical 
predictions recently given by Castelnovo et al. In the parameter range where the theory can be 
applied we find almost perfect quantitative agreement. 
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I. INTRODUCTION 



The term polyelectrolyte is employed for the wide field 
of macromolecules which contain dissociable subunits 
(see, e.g., [[j], ^ and references therein). With respect 
to different dissociation behavior one can distinguish be- 
tween strong and weak polyelectrolytes [Q or between 
quenched and annealed ones S. So-called strong poly- 
electrolytes, polysalts as e.g. Na-polystyrene sulfonate, 
dissociate completely in the total pH range accessible by 
experiment. The total charge as well as its specific dis- 
tribution along the chain is solely imposed by chemistry, 
i.e., by polymer synthesis. That is why such polyelec- 
trolytes are also called quenched. On the other hand, 
weak polyelectrolytes represented by polyacids and poly- 
bases dissociate only in a rather limited pH range. The 
total charge of the chain is not fixed but it can be tuned 
by changing the pH of the solution. Because of dissocia- 
tion and recombination of ion pairs along the chain one 
expects spatial and/or temporal fluctuations in the local 
degree of dissociation. Such titrating polyelectrolytes ex- 
hibit an annealed inhomogeneous charge distribution. A 
pronounced charge accumulation appears at chain ends 
because there are fewer neighbors for the charges to inter- 
act with and the penalty in energy is therefore reduced. 
Although, at the level of scaling laws describing the sta- 
tistical properties of polymer chains, the local charge dis- 
tribution has only a weak effect on numerical pre-factors 
H the extra degree of freedom for the charges leads to 
new and non trivial features. The charge inhomogene- 
ity can have a strong impact on processes dominated by 
end-effects, such as the self-assembly of weakly charged 
linear micelles B and adsorption on charged surfaces 
||. For end-grafted weak polyelectrolytes, a rather un- 
usual regime has been obtained where the chain stretch- 
ing (brush thickness) depends non-monotonously on salt 
concentration || and grafting density |7| . This is mainly 



due to the fact that the net charge of a chain as well as 
its distribution along the chain is not fixed but depends 
on its local environment. 

The dissociation of a low molecular acid (HA) in an 
aqueous medium is given by the equilibrium reaction 

HA + H 2 ^ H 3 0+ + A - , 

or more simply by 

HA ^ H+ + A". 
The law of mass action yields the equilibrium constant 

[H+] [A-] 



K„ 



[HA] 



(1) 



where [A - ], [HA], [H + ] are the (monomolar) concentra- 
tions of dissociated and undissociated acid and dissoci- 
ated hydrogen, respectively. Using the standard notation 
pH = -logio [H+], pK a = -logio [K a ] and defining the 
degree of dissociation by 



/ 



[A- 



[HA] + [A- 



(2) 



eq. ([!]) gives the well-known relation between the pH of 
the solution and degree of dissociation / of a simple acid 



pH = pK a + log- 



io 



1-/ 



(3) 



The dissociation behavior of polyacids can be described 
in a similar way, but the resulting pK a value is now an 
apparent one (in the physico-chemical literature denoted 
by pK app ) ||. In contrast to low-molecular- weight acids, 
the charged groups of polyacids are linked together along 
the chain. Therefore, the dissociation of one acid group 
is correlated in a complex way to the position and the 
number of other charged groups of the chain resulting 
in a masking of the intrinsic pK® of a (polyelectrolyte) 
monomer. The corresponding relation can be written 
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P H = pK° + log- 
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dF el 
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with F e i being the electrostatic free energy of the poly- 
electrolyte chain and iV is the chain length. Introducing 
the chemical potential 



eq. 



= fcsTlogi 



can be written 



1-/ 



1 dF a 



n df 



pH = P K° a 



(5) 



(6) 



Clearly, the chemical potential has two contributions: (i) 
an entropic one related to the mixing of charged and 
non-charged groups along the chain and (ii) an electro- 
static one related to the interaction with charged groups 
forming the local charge environment of an ionizable site. 
For good and 9 solvents, the electrostatic contribution 
A t ei(/) = N^ 1 dF e i/df is an increasing monotonic func- 
tion of /. For poor solvents, Raphael and Joanny |J 
found a non-monotonic variation of ^ e ;(/), and thus of 
fj-(f), with / which results in a first order phase transi- 
tion from a collapsed weakly charged conformation to an 
extended strongly charged state. 

The particular behavior of weak polyelectrolytes has 
attracted considerable interest in experimental Efl 1 1 . 
ll, theoretical || |J, |l§ and simulation fTj, [17|, 
18, |l9[ |2^] studies. Because a usual experimental ap- 
proach to characterizing weak polyelectrolytes is to per- 
form titration experiments much effort has been done to 
understand the titration curves, i.e., the dependence of 
the degree of ionization (or the degree of neutralization in 
the titration experiment) on the pH of the solution. The 
inhomogeneous charge distribution has been first studied 
by numerical simulation |l4|, |l9) . For not too large mean 
charge densities (/) and screening lengths Xd, a quite 
good agreement between simulation data and the predic- 
tions of a linearized density functional approach has been 
obtained in the case of rod- like polyelectrolytes jl| fli) . 
Recently a generalization of the theory to the case of 
flexible chains has been given by Castelnovo et al. Jl6| . 
Qualitatively they found the result obtained in our pre- 
vious simulations: a charge accumulation at chain ends. 
However, a quantitative comparison between theory and 
simulation data was not possible because the chains con- 
sidered by Berghold at al. |l4| are not in the asymptotic 
regime where the approach used in theory can be applied. 
A similar end-effect has been demonstrated in simula- 
tions of quenched strongly charged polyelectrolytes with 
explicit counterions: around the ends the counterion dis- 
tribution is significantly different from that around the 
inner part of the chain However, although the re- 

sulting effective charge looks quite similar to that of an 
annealed polyelectrolyte a quantitative comparison with 
the corresponding theoretical predictions is not possible 
due to some important differences between the two sys- 
tems. So the theoretical predictions given by Castelnovo 
et al. are still waiting for a quantitative verification. 

In the present paper we study the non-uniform charge 
distribution on annealed weakly charged polyelectrolytes 



in a solvent by (semi-)grand canonical Monte Carlo sim- 
ulation. To be close to the theoretical model we consider 
chains where neighboring monomers are bound by har- 
monic springs. In addition, charged monomers interact 
with a Debye-Huckel potential, the screening length of 
which is tuned over the range where the theory is valid. 
In this case we observe not only qualitative agreement 
between simulation data and theory, but indeed a quite 
good quantitative one. Although the theoretical single 
chain problem and its solution in a restricted parameter 
range may seem rather academic it is a first step towards 
a more precise understanding of unusual collective prop- 
erties of, e.g., annealed polyelectrolyte brushes and stars. 

The outline of the paper is as follows. In Section || 
we present the main theoretical predictions. The model 
and the method we use in the simulation are described 



in Section III. In Section IV we discuss the results and 



compare theory and simulation. Finally, a brief summary 
and our conclusions can be found in Section |v|. 



II. THEORY 

The free energy of an annealed, fully-stretched poly- 
electrolyte chain [N monomers of size b) in a salty solu- 
tion can be written |E3| 



F 



N/2 



N/2 



dJ/(a) [log (/(*)) - 1] - M/W 



Ab f N/2 .expHs-s'l&M^) 



+ ds'f(s)f(s') 



N/2 



\s-s'\b 



(7) 



where f(s) is the local charge distribution along the 
chain. The Bjerrum length, which sets the strength of 
electrostatic interactions, is given by A^ = e 2 /A-KE ekBT 
with e being the elementary charge and e is the di- 
electric constant of the medium. Assuming a mono- 
valent low- molecular electrolyte of concentration c s , De- 
bye screening length Ad and Bjerrum length arc related 
by \ D = (8ttX b c s )' 1/2 . The first term of eq. (@) is the 
entropy of an unidimensional ideal gas, the second term 
fixes the charge on the chain by the chemical potential 
fi and the third term represents the electrostatic inter- 
action between the charges on the chain. Note that one 
can justify the use of a Debye-Huckel potential only if 
the average charge density (/) / b fulfills certain condi- 
tions: (i) (/) has to be small to make non-linear effects 
unimportant, i.e., 



</>^ = (/>« < 1, 



(8) 



where we introduced the dimensionless parameter u = 
As/fr, and (ii) (/) has to be large to make sure that a 
sufficiently large number of charges interact simultane- 
ously, i.e., 



(/) ^ > 1. 



(9) 
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Note that (/) Xs/b is the so-called Manning parameter 
for the condensation of counterions on a partially charged 
rigid rod. 

Minimizing eq. (R) the equilibrium charge density dis- 
tribution on an annealed, fully stretched polyelectrolyte 
chain, up to first order in (/) Xs/b, was found to be 



(f) x f 



(10) 



N 
2 



N 

2 ~ S J X D 



2^ 
Nb 



with the exponential integral Ei(x) = dtt -1 exp(— t) 
p3| . Equation (10) gives a charge accumulation at the 
ends of the rod within a region of the order of the screen- 
ing length. For a wide parameter range, the theoretical 
prediction was quantitatively confirmed by Monte Carlo 
simulations pi . 

In order to study the charge distribution on a weakly 
charged, flexible polyelectrolyte chain Castelnovo et al. 
fL6f generalized eq. (Q) by including the entropy of a 
freely jointed chain. In the 9 solvent case the free energy 
can be written 



F[r (s) 



N/2 
■N/2 

3 

262 



ds\f(s) [log (/(a)) - 1] - 



dr (s)\ 2 , X B rm 



ds 



ds'f(s)f(s') 



-N/2 



exp[-|r (s) - r (s')\/X 1 
|r (*)-ro(8')l 



(11) 



(36) - (39)), for a weakly charged chain of electrostatic 
blobs the conditions for applying the Debye-Hiickel ap- 
proximation, which correspond to eqs. (||, ||) discussed 
above, become 



and 



As£o 



b 2 



A.p£o 

b 2 



< 1, 



> 1. 



(14) 



(15) 



Thus, the Manning parameter reads now (/) Xb £,o/b 2 ~ 

1/3 

((/) u 2 ) 1 . Note that the condit ions imposed by apply- 
ing the chain under tension model, i.e., (i) a fully flexible 
chain with £ 3> b and (ii) a aligned blob chain with 
Xn £o, are even stronger than those given in eqs. ( p^ , 
|l5| ) . Specifically, assuming that 



b <C Co < Ad < L 



V 



(16) 



where L p is the persistence length (for inherently flexible 
chains due to electrostatic interaction), the solution of 
f(s) up to first order in (J) Xb t;o/b 2 was found to be 
very similar to eq. (jLCfl) p6| 



M 



Xr 



(17) 



N 
2 



Xi 



N 



A 



D. 



Nb 



where ro(s) denotes the so-called classical path of the 
polymer (the most probable conformation) which follows 
from the "chain under tension model" |2j] . An important 
length scale which comes now into play is the electrostatic 
blob size £■ From scaling arguments, i.e., by equating the 
(unscreened) electrostatic interaction of the (quenched) 
charges inside one blob with ksT, one gets 



^scaling 



(upf 3 



(12) 



One blob contains g — (£scaiing/&) 2 
of TV monomers consists of 



monomers and a chain 



n b = N/g = N(uf 2 ) 2/3 



(13) 



electrostatic blobs. In order to reach the asymptotic 
regime one has to ensure that nj, 3> 1. For a chain of finite 
length, however, the situation becomes even more com- 
plex. The electrostatic blob size varies along the chain 
resulting in a trumpet-like conformation with the maxi- 
mum blob size at the chain ends. Up to a factor A the 
minimum blob size £o occurring in the middle of the chain 
is equal to ^scaling 10 ■ Later on we will use A as an ad- 
justable parameter for comparison with simulation data. 
Using the rescaling relations given in reference J16[ (eqs. 



with Xb,Xd being Bjerrum and screening lengths ex- 
pressed in terms of contour length 



Ar 



3£o 



A_D — —r- AD- 




(18) 



(19) 



The relation between persistence length and screening 
length is a controversial problem which has been dis- 
cussed for a couple of decades since the early work of 
Odijk, Skolnick and Fixman [^5[ ||(|. It is beyond the 
scope of the paper to go into details of that problem. 
Note that for fully flexible chains, the most recent re- 
sults confirm the prediction for a chain of electro- 
static blobs originally made by Khokhlov and Khacha- 
turian Ip8l 



Sscaling 



4 ^scaling 



(20) 



Below we will use this relation to estimate the persistence 
length of the weakly charged chains considered in the 
simulations (see Table 1). 



4 



TABLE I: Systems studied and related scaling quantities 
(N = 1000, u = X B /b = 0.9). 
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16 
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1.5 


0.040 


8.9 


13 


0.32 


64 
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3.7 










256 


28.8 


16.7 










16 
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1.8 


0.083 


5.4 


34 


0.40 


64 


11.9 


6.3 










256 


47.4 


25.0 










16 


3.9 


2.2 


0.125 


4.1 


58 


0.46 


64 


15.6 


8.3 










256 


62.4 


30.9 



III. SIMULATION MODEL AND METHOD 

In the simulation the polyelectrolyte is represented 
by a freely jointed bead-spring chain, the realization of 
which we chose as close as possible to the model used in 
theory. Along the chain the N monomers are connected 
to their neighbors by a harmonic bond potential 



t^bond = t,^bT ^ 
1 n=l 



N-l , 

( r n+l 



Ir 



(21) 



with r„ being the position of bead n and bo is the (bare) 
average bond length, henceforth set equal to one. The 
thermal energy is fc^T. For convenience, we use ksT = 1. 
All N c (negatively) charged monomers interact via the 
Debye-Hiickel potential 



■ E 

n^m—l 



■ exp 



' nm 

X7 



(22) 



where the Debye screening length Xd is an input param- 
eter. The exponentially decaying interaction enables the 
introduction of a cutoff which we chose as A c = 5 . 
In the case of weakly charged polyelectrolytes in a 6 sol- 
vent, we study in the paper, average bond length as well 
as bond length distribution are only slightly effected by 
the electrostatic repulsion between charged monomers. 
Thus, the average bond length remains b « bo = 1. For 
water at room temperature, the Bjerrum length which 
gives the strength of the Coulomb interaction is about 
7.14 A. To avoid problems with counterion condensa- 
tion, and to ensure that we work in a parameter range 
where the theory discussed above can be applied, we set 
the length scale by u = Xs/b = 0.9. Then the Man- 
ning parameter of a chain of electrostatic blobs (see eq. 

@ ) obeys (/) \ B io/b 2 ~ ((/) u 2 ) 1/3 < 1 by definition. 
With this setting of the length scale one has b « 8 A. 
Hence, the polyelectrolyte chain is modeled on a coarse 
grained level where one bead corresponds to a few chem- 
ical monomers. In the simulations reported here we con- 
sider a chain which consists of N = 1000 beads N c of 



which carry an elementary charge. Because we study an 
annealed chain in a (semi-)grand canonical ensemble, N c 
is not constant but only its average value (N c ) is fixed 
by the chemical potential \x. The chemical potential is 
chosen to result with an average degree of dissociation 
(/> = {N c ) /N ranging from 1/25 to 1/8. The Debye 
screening length A^> is varied over a range from 16 to 256 
average bond lengths. With the setting of length scale 
introduced above such screening lengths correspond to 
salt concentrations from 10 _3 mol/L to 10 _5 mol/L. Ta- 
ble 1 shows the systems studied by simulations together 
with the scaling quantities being of interest with respect 
to the conditions given in eq. (1^). One can see that the 
polyelectrolytes considered in the simulations at least ful- 
fill the inequality relations which ensure that the chains 
reach the asymptotic regime assumed in theory. The con- 
ditions ri), 3> 1 and (/) \b /b <C 1 discussed in Section 2 
are also reasonably fulfilled. 

Equilibrium properties of the polyions are studied by 
standard Metropolis Monte Carlo (MC) |^9| simulation. 
In order to guarantee the equilibration of both long and 
short length scales, bond angles and bond lengths we 
combine two different configurational MC moves: (i) A 
pivot move where a monomer is chosen at random and the 
subsequent part of the chain is randomly rotated around 
that monomer, and (ii) a local displacement move. The 
pivot algorithm was shown to be a highly efficient way to 
sample the phase space in a single chain model with fixed 
bond lengths M. In the case of our freely jointed 
bead-spring model we found that the most efficient sim- 
ulation route is a 1:1 mixing of the two MC moves. The 
correlation time of the mean square end-to-end distance 
has been checked to be of the order of a few tens of MC 
steps. Thus, in any case we have correlation times below 
0.1 Monte Carlo steps per monomer (MCM). 
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FIG. 1: Chemical potential /i versus average degree of dis- 
soziation {/) for a completely extended rigid chain at dif- 
ferent Bjerrum lengths: As =0 (solid line), As = 0.56 
(dashed), As = 1-06 (dot-dashed). Simulation results with 
N = 100, X D = 106. 



In order to simulate annealed polyelectrolytes the MC 
simulation is performed in a semi-grand canonical en- 
semble where the chain is in contact with a reservoir of 
charges of fixed chemical potential /i. Additionally to 
the configurational MC moves introduced above the al- 
gorithm is completed with a charge move by which the 
charge state of a randomly chosen monomer is switched. 
Again a 1:1:1 combination of the three different MC 
moves was found to be the most efficient choice. The 
energy change of a complete MC move reads 

AE = AE C ± fi, (23) 

where AE C is the change in configurational energy due to 
^bond and C/dh- The plus sign is used when the monomer 
is to be neutralized (protonated) and the minus sign when 
it is to be charged (deprotonated) . Figure |l| shows the 
dependence between \i and average degree of dissociation 
(/) at different strength of Coulomb interaction, obtained 
by simulating a fully extended (rigid) polyelectrolyte. At 
vanishing interaction we have the ideal curve of isolated 
monomers. According to eq. (0) /i/ksT equals pH— pif°. 
Thus, Figure |l] represents the titration curves of the spe- 
cific model. It might seem rather academic to study the 
behavior at different Bjerrum lengths. However, the im- 
portant dimensionless parameter giving the strength of 
Coulomb interaction is u = Xs/b. Hence, varying the 
average distance between ionizable groups 6, which can 
be easily done during polyelectrolyte synthesis, one can 
succeed with a similar tuning of the strength of interac- 
tion as by varying A^. 

Exploring different starting configurations, specifically 
a random one and the completely stretched one, the equi- 
libration time of the chain was estimated to be r equ ss 500 
MCM. To ensure sufficient relaxation, first the simula- 
tion is run for a time of at least 160 r equ which are about 
80 • 10 6 MC steps. After reaching equilibrium, ensemble- 
averaged chain properties are taken as averages over at 
least 1.6T0 4 MCM, which corresponds to at least 1.6T0 5 
renewal times. Note that we consider one MC step to be 
completed when the three partial steps (Pivot move, lo- 
cal move, charge move) have been accepted. In this way 

TABLE II: Conformational properties of the systems studied: 
end-to-end distance R, radius of gyration R g and shape factor 
R 2 /R 2 g {N = 1000, u = \ B /b = 0.9). 
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X D 
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16 


62.8±0.3 


23.6±0.1 


7.1±0.1 


0.040 


64 


75.2±0.3 


27.3±0.1 


7.6±0.1 




256 


77.4±0.3 


28.0±0.1 


7.6±0.1 




16 


98.0±0.4 


35.5±0.1 


7.6±0.1 


0.083 


64 


132.9±0.4 


45.5±0.1 


8.5±0.1 




256 


143.5±0.3 


48.5±0.1 


8.8±0.1 




16 


129.3±0.4 


46.1±0.1 


7.9±0.1 


0.125 


64 


184.5±0.4 


61.5±0.1 


9.0±0.1 




256 


206.2±0.4 


67.4±0.1 


9.4±0.1 




FIG. 2: Simulation snapshots of a partially charged chain 
(iV = 1000, regularly quentched distribution with / = 1/8, 9 
solvent) at a) Ad = b and b) Xd = 500 6. 

the results on relevant quantities were found to be repro- 
ducible within a few percent. 

Note a significant dependence of the acceptance rate 
on (/) and A_d. We estimated it to range from about 
60% for (/) = 0.040, X D = 16 6 down to about 20% for 
(/) = 0.125, Xd = 256 6. Due to this difference CPU 
times between 16 hours and 9.5 days were necessary for 
performing simulation runs on Compaq Alpha machines 
with EV67/667MHz processors. 



IV. SIMULATION RESULTS AND 
COMPARISON WITH THEORY 

Table || shows conformational properties of the chains 
studied by simulations. As expected the stretching of 
the chains is growing with increasing charge fraction (/) 
and/or screening length Xd- The largest shape factor 
R 2 /R 2 g we obtain is 9.4. Remember that it ranges from 
6 for a Gaussian chain up to 12 in the case of a rigid rod. 
For the largest screening length Xd = 256 6, we have al- 



6 



ways Xd > R, i.e., the electrostatic interaction is almost 
unscreened. Note that the effect of fluctuations in the 
charge distribution on large-scale mean statistical prop- 
erties of free polyelectrolyte chains in dilute solutions is 
weak. In Figure ^ two typical simulation snapshots of 
a partially charged chain are shown where unite charges 
are regularly distributed with / = 1/8. The qualitative 
difference between the two snapshots, i.e. the polyelec- 
trolyte effect, is not effected by charge annealing. In the 
high salt case a), the Coulomb interaction is almost com- 
pletely screened (Xd — b) and the configuration looks 
quite similar to a swollen coil in a good solvent. On 
the other hand, in the weak screening case b) we have 
an almost linear chain of blobs. At short length scales, 
i.e. inside the blobs, there are coiled regions. At large 
length scales the blobs appear to build a strongly elon- 
gated chain the transverse extension of which is consider- 
ably smaller than the longitudinal one. An appropriate 
quantity that describes quantitatively the structure of 
the chain at all length scales is the single chain structure 
factor or form factor. Here we calculate the spherically 
averaged structure factor 



S(q) 




N 

E 

n=l 



exp(zq • r„ 



(24) 



From the theory of uncharged polymers we know that 
the structure factor scales as S(q) ~ q~ x l v in the range 
2tt/R < q < 2n/b with v being the universal exponent 
for the mean extension of the chain R ~ N v . The ideal 
chain and good solvent chain values of v are 1/2 and 
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FIG. 3: Spherically averaged structure factor for partially 
charged chains (N — 1000, annealed charge distribution). 
Simulation results at varying degree of charging and screening 
length: (/) = 0.040, \ D = 16 6 (circles); (/) = 0.083, X D = 
64b (squares); {/) = 0.125, X D = 256 6 (triangles). Addition- 
ally the result of an ideal chain is plotted (diamonds). The 
thin lines indicate asymptotic scaling laws. 



0.588 (w 3/5), respectively. In Figure | S(q) is plotted 
for three different systems with annealed charge distribu- 
tion studied by simulation: (i) At the minimum extension 
of the chain, (ii) at a mean one and (iii) at the maximum 
one (see Table 0). For large degree of charging and not 
too short screening lengths linear scaling with v « 1 is 
reached at large length scales. On short scales we have 
an almost ideal random coil behavior. This is exactly 
the structure implied by the theory. Thus, in contrast 
to the previous simulations by Berghold et al. Jl4| now 
we expect a rather good agreement between simulation 
data and theory. In Figure we show the corresponding 
charge distributions together with the theoretical predic- 
tions following from eq. ( |l7| ) . We remember that the the- 
ory contains a free parameter A which sets the relation 
between the scaling theory blob size £ SC aiing and blob size 
in the middle of a finite chain £o = A. Scaling- We fit A to 
1.5- 



1.4- 



1.3 




1.3 



0.9 



a) 



I I I I I I I I I I — \ 




b) 




i 1 1 1 1 i 1 1 1 1 i 




c) 





-0.5 0.5 

nl{Nl2) 



FIG. 4: Equilibrium charge distribution on annealed par- 
tially charged polyelectrolytes (N — 1000). Simulation results 
(symbols) and theoretical predictions (lines) at varying degree 
of charging and screening length: a) (/) = 0.040, Ad = 16 6, 
b) (/) = 0.083, A D = 64 6, c) (/) = 0.125, X D = 256 6. 
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TABLE III: Fitting of charge distribution /(s) with the 
charge density in the middle of the chain /(0). Values of 
parameter A. 
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0.040 
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16 


0.083 


1.8 




0.125 
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0.040 


1.3 


64 


0.083 


1.2 




0.125 


1.2 




0.040 


0.89 


256 


0.083 


0.80 




0.125 


0.76 



obtain best agreement of the charge density at the middle 
of the chain where it forms an exten ded plateau. The re- 
sulting values of A are given in Table 



III 



In doing so, the 
agreement between simulation data and theoretical pre- 
diction is indeed almost perfect for the case where all the 
three inequalities given in eq. (16) are fulfilled (case b) in 
Figure Q). For decreasing charging (case a) in Figure^), 
the size of electrostatic blobs grows and becomes of the 
same order of magnitude as screening length (see Table 
||) . Then the interaction is too strongly screened to align 
the blob chain. At (/} = 0.040, \ D = 16 6 the structure 
factor clearly shows (see Figure ||) that the chain docs 
not reach the asymptotic regime with v « 1 assumed in 
the theoretical model. Because of the stronger coiling a 
longer part of the chain is packed within the range of a 
screening length which sets the spatial scale of the charge 
inhomogeneity. Thus, in agreement with simulation data 
we expect the charge to be accumulated in a larger chain 
section than predicted by theory. For strong charging 
(case c) in Figure^), the expansion parameter (/) Ad/6 
becomes too large to justify a first order perturbational 
treatment used in the theoretical approach discussed in 
Section O. At this point we have to note that the fitting 
procedure of A described above is actually more complex 
than simply finding the right pre-factor of Scaling in the 
case of finite chain length. We remember that due to the 
restrictions of a first order theory the normalization of 
the charge distribution f(s) given by eqs. (|l(], [l7]) is cor- 
rect in the limit N ^> Ab/6 or N 3> Ad/6, respectively. 
In reference jl4j it has been shown that the first order 
corrections due to finite size effects give an additive con- 
stant in f(s). Thus, the fitting of the free parameter A 
corrects not only for the unknown pre-factor of Scaling 
but also for the higher order terms neglected in theory. 
However, doing so we obtain an almost perfect agreement 
between simulation data and theoretical predictions for 
the charge accumulation at chain ends within a length 



of the order of the Debye-Huckel screening length. From 
eqs. (|l7|), ( [l8|) we know that, in first order theory, the 
amplitude of charge inhomogeneity is proportional to the 
rescaled strength of interaction, i.e., among others it is 
proportional to £o- On the other hand, Table III shows 
that the relation between £o and Scaling is monotonously 
reduced with growing Ad and/or (/}. Thus, the larger 
screening length and degree of charging, the stronger be- 
comes the correction of the (straightforward) first order 
result. Actually this tendency is in complete agreement 
with the behavior we obtained for rigid rods without in- 
troducing any fitting parameter With growing Ad 
and / the first order theory overestimates the inhomo- 
geneity, but this error is compensated (at least partially) 
by the first correction term. 



V. CONCLUSION 

In this paper we reported (semi-)grand canonical 
Monte Carlo simulations of annealed weakly charged 
polyclcctrolytcs. In a fairly wide parameter range, i.e., 
not too large mean charge densities (/) < 6/Ad and 
screening lengths in the range Scaling < Ad < Nb we 
find a quite good quantitative agreement between sim- 
ulation data and the results of the linearized theory re- 
cently proposed by Castelnovo et al. (T(|. It would be 
interesting to have experimental evidence for the accu- 
mulation of charge at chain ends we obtain within a re- 
gion of the order of the screening length. In order to 
being able to compare simulation data with theory we 
have restricted the simulations to weakly charged poly- 
electrolytes ((/) < 0.125). However, from previous sim- 
ulations |l4j it is known that the maximum of the over- 
charging at chain ends appears close to (/) rs 0.5 where 
the charge density at the ends can become about 50% 
higher than its mean value (/}. In this case it should be 
possible to see experimentally the effect of the charge in- 
homogeneity on processes dominated by end-effects such 
as, e.g., self-assembly of weakly charged micelles jl and 
adsorption on charged surfaces ||. But, such a charge 
density is clearly outside the range where the first order 
theory is valid. Nonetheless the restricted theory, now 
for the first time proved by simulation data, is a step to- 
wards a more precise understanding of the more complex 
systems mentioned above. 
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